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ABSTRACT 

We propose that ultra-high-energy (UHE) cosmic rays (CRs) above 10^® eV are produced in rel¬ 
ativistic jets of powerful active galactic nuclei via an original mechanism, which we dub “espresso” 
acceleration: “seed” galactic CRs with energies < 10^^ eV that penetrate the jet sideways receive a 
“one-shot” boost of a factor of ~ E^ in energy, where T is the Lorentz factor of the relativistic flow. 
For typical jet parameters, a few percent of the CRs in the host galaxy can undergo this process, 
and powerful blazars with T > 30 may accelerate UHECRs up to more than 10^° eV. The chemical 
composition of espresso-accelerated UHECRs is determined by that at the Galactic CR knee and is 
expected to be proton-dominated at 10^® eV and increasingly heavy at higher energies, in agreement 
with recent observations made at the Pierre Auger Observatory. 

Subject headings: Acceleration of particles — cosmic rays — galaxies: active 


1. THE ORIGIN OF ULTRA-HIGH-ENERGY (UHE) 
COSMIC RAYS (CRS) 

The mechanism responsible for the acceleration of 
highest-energy CRs detected at Earth is one of the 
most prominent open questions in astrophysics. Possi¬ 
ble sources of particles as energetic as I0^° eV are lim¬ 
ited to gamma-ray bursts (GRBs) (e.g., Vietri 1995; 
Waxman 1995), newly born millisecond pulsars (e.g., 
Blasi et al. 2000; Fang et al. 2012), and active galactic 
nuclei (AGNs) (e.g., Ostrowski 2000; Murase et al. 2012), 
which can —in principle— generate the electric poten¬ 
tial necessary to produce UHECRs (Hillas 1984; Dermer 
2007). Recent improved measurements of extensive at¬ 
mospheric air showers induced by UHECRs are providing 
unprecedented information on their chemical composi¬ 
tion. The Pierre Auger Observatory measured a compo¬ 
sition compatible with pure protons around 10^® eV, and 
increasingly heavy at higher energies (Aab et al. 2014a). 
Such results are not inconsistent with Telescope Array’s 
observations (Abbasi et al. 2015), which favor a lighter 
chemical composition, especially if interpreted in light of 
the LHC data (Pierog 2013). The evidence for nuclei 
with charge Z > 1 in the most energetic UHECRs may 
be of paramount importance for unraveling their origin. 

In this Letter, we exploit some basic properties of par¬ 
ticle acceleration in relativistic flows to outline a simple 
but general mechanism for the production of UHECRs in 
powerful AGNs. We call such a process the “espresso” 
mechanism, since a very fast flow (the AGN jet) pro¬ 
vides a one-shot acceleration of energetic seeds (galactic 
CRs). The implications of such a scenario for CR spec¬ 
trum and chemical composition at energies E > 10^® eV 
are discussed and compared with observations. 

2. AGNS AS POTENTIAL SOURCES 

We briefly review the requirements for candidate 
UHECR sources (energetics, luminosity, and confine¬ 
ment) and check that AGNs satisfy all of them, if the 
highest-energy CRs are heavy nuclei. 

The energy generation rate necessary to sus- 
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tain the UHECR flux above lEeV is Ri5.4 x 
10^®ergMpc“®yr“^(e.g., Katz et al. 2013). Typ¬ 
ical bolometric AGN luminosities range between 
Lboi «10^®ergs“^ for Seyfert galaxies and radio-quiet 
quasars, to Lboi ^ lO^^ergs"^ for radio-loud quasars 
(e.g., Woo & Urry 2002); since the typical density of ac¬ 
tive galaxies is « 10“''’Mpc“®, 10”20% of which is radio- 
loud (Jiang et al. 2007), the total AGN energy genera¬ 
tion rate (in photons) is ss5xlO"‘®ergMpc“®yr“^. AGNs 
would meet the energy requirement if they emitted a frac¬ 
tion > 10“® of their bolometric luminosity in UHECRs. 

If the flow is expanding, the rapid acceleration up to 
an energy iKmax may also require a minimum source lu¬ 
minosity (Waxman 2004): 
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( 1 ) 


where P and jdc are the flow Lorentz factor and speed, 
and ^26 is the particle charge in units of the iron nucleus 
charge. For powerful AGNs, Lboi ^ A,, provided that 
UHEGRs with Umax ~ 10^° eV are heavy nuclei and that 
r < 50. However, Lboi represents a lower limit on the 
actual source power, and the luminosity in CRs can be 
much larger than that in photons: for instance, Lcr ~ 
lO^Urgs-i » Lboi < 10®®ergs-i in Tycho’s supernova 
remnant (SNR) (Morlino & Caprioli 2012). In blazars 
the jet power is inferred to be a factor of 10-100 larger 
than Lboi (e.g., Ghisellini et al. 2009). 

The constraint in Equation 1 is similar (though not 
equivalent) to the so-called Hillas criterion (Hillas 1984), 

A;.G-Rkpc > io20eV’ 

which expresses the maximum energy achievable in a sys¬ 
tem with magnetic field B (in jjtG) and size R (in kpc). 
For typical AGN jets, U^cAkpc ^ 1, allowing confine¬ 
ment and acceleration of nuclei up to U > 10^° eV. 

These considerations support the hypothesis that 
AGNs can accelerate EeV protons and iron nuclei up 
to the highest observed energies even for rather small 
acceleration efficiencies (> 10“^ of the jet power). 
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3. ACCELERATION IN RELATIVISTIC FLOWS 

Non-relativistic supersonic flows that dissipate their 
energy into shocks can efficiently energize CRs via diffu¬ 
sive shock acceleration (e.g., Caprioli & Spitkovsky 
2014a); in relativistic flows, instead, diffusive 
acceleration may be quite suppressed (see, e.g., 
Sironi & Spitkovsky 2011), and different mechanisms 
are needed to produce energetic particles. 

Let us consider a relativistic flow with Lorentz factor F 
and velocity /3c x in the laboratory frame, and a particle 
with initial energy and momentum 

Pi ~ /^?,0), (3) 

where p = Px/\p\', once in the flow, its energy in the flow 
frame is 


E! = TiEi - /3pi^,) = TEiil - /3pi). (4) 


If the particle gyrates around the comoving magnetic 
field B' before leaving the flow, its final energy and flight 
direction can be written as Ef = E'^ and p( = Pf ^/E[^ 
which in the laboratory frame become 

Ei=T^E,(l-Pp{){l + Pp[), (5) 

If Pf = pi^ the particle energy is unchanged, but typ¬ 
ically Pf 7 ^ Pi, which implies E{ ~ F^ifi, similar to a 
Compton scattering against a relativistic magnetic wall. 
This phenomenon, which is independent of where par¬ 
ticles enter/leave the flow, is well known for relativistic 
shocks: the energy gain is ^ 2F^ in the first upstream- 
downstream-upstream cycle (pi ~ — l,Pf — I), but ^ 2 
in the following ones because particles are re-caught by 
the shock with p > 1 —1/F ^ pf (Achterberg et al. 2001). 

Let us consider a particle with yin ^ F entering the 
flow with Pi = 0, corresponding to p( = —/3, and assume 
B' = —B'z. In the flow frame, the particle performs a 
Larmor gyration with frequency ^l' = eB'/ me) and 

in turn (Equation 5): 


Ei 


1-/3^ cos n't' + ^ sin n't' 


( 6 ) 


The total energy gain depends on the phase ip' = n't' 
(mod 27r) when the particle leaves the flow: it is maxi¬ 
mum (2F^) for (pf = 7r/2, and ~ F^ for p'f G [7r/2, 37r/2]. 
A boost of ~ F^ in the laboratory requires the particle 
to stay in the flow for Tacc ^ r7r/(2n'), during which it 
travels a distance 

i^acc TaccC « 4kpC^^^^|^-, (7) 

with 7 f = 7 inF^ ^ 5 X 10® the maximum UHECR Lorentz 
factor. In reality, since relativistic flows diverge and the 
magnetic field drops {B' few G a:";) in blazar jets, 
e.g., O’Sullivan & Gabuzda 2009), particles eventually 
escape, either because they reach the termination shock 
or because the condition in Equation 2 is violated. Also, 
the expected radial dependence of the toroidal magnetic 
field induces an axial VS—drift toward the flow head. 

In realistic flows, we can assume that the orbit is gener¬ 
ally truncated with a random phase p'f , which leads to an 



Figure 1. Schematic trajectory of a galactic CR reaccelerated by 
a relativistic jet (not in scale). The total acceleration due to the 
motional electric field does not depend on the exact trajectory: a 
rotation ^ 7r/2 around the jet magnetic field is sufficient to achieve 
a boost (Equation 6). 


average energy gain (Ef/Ei)^'^ = F® (Equation 6). Fig¬ 
ure 1 shows the sketch of a possible particle trajectory in 
a conical (expanding) jet. Exact trajectories in realistic 
velocity/magnetic profiles of AGN jets will be presented 
in a forthcoming publication, but we anticipate that ~ F® 
energy gains are indeed common in astrophysical jets. 

UHECR acceleration via such a one-shot (espresso) 
mechanism thus requires either ultra-relativistic flows 
with F > 10®, or moderate Lorentz factors and pre¬ 
accelerated particles. We now consider the case of AGN 
jets reaccelerating energetic CR seeds. 


3.1. “Seeds” 

Galactic CRs accelerated in SNRs (Morlino & Caprioli 
2012; Ackermann et al. 2013) represent natural seed can¬ 
didates. The maximum energy Amax ~ few Z PeV (the 
CR “knee”) is achieved before the SNR enters the Se¬ 
dov stage (t ss ts), when the shock velocity Us starts to 
decrease because of the inertia of the swept-up material 
(Blasi et al. 2007). E^^x can be estimated by equating 
ts and the acceleration time face which scale as 


tacc OC 


Ej 

-^max 

W’ 






( 8 ) 


where Esn and Mej are the SN ejecta kinetic energy and 
mass, and p and B are the circumstellar density and mag¬ 
netic field; we have also used Bohm diffusion to derive 
tacc (see Caprioli & Spitkovsky 2014b,c, for a justifica¬ 
tion of this assumption based on ab-initio simulations). 
Finally, we have 

E^ay. OC ZBp-fe® OC Zp^/'^^/T~, (9) 

where we also assumed equipartition between thermal 
and magnetic pressures (as in the Milky Way), i.e., oc 
pT, with a typical temperature T of the order of the 
virial temperature Tvir. Since the dependence on p is very 
weak, and since Tvir is proportional to the total galactic 
mass, which does not differ greatly from galaxy to galaxy, 
Bmax is expected to be roughly the same for any SNRs 
expanding in the interstellar medium. For core-collapse 
SNe, instead, Bmax is achieved while the shock is still 
propagating in the wind launched in the pre-SN stages 
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(Bell et al. 2013; Cardillo et al. 2015), which should be 
independent of the properties of the galaxy. We conclude 
that the maximum energy of CRs accelerated in SNRs 
should be rather universal and correspond to the CR 
knee in the Milky Way. 

3.2. “Steam” 

Let us consider a galaxy that hosts a powerful active 
nucleus launching a relativistic jet with E Ri 30 and an 
opening angle of Ad Ri 2°, which propagates for Hj ^ 
several kpc through diffuse galactic CRs. The typical 
gyroradius of knee nuclei (~ Ipc Ep^y/ZB~ q) is much 
smaller than the transverse size of the jet (i?j « HjAd), 
but larger than the jet boundary layer, whose thickness 
is determined by small-scale plasma processes. 

We now estimate the fraction of galactic CRs that 
can percolate through the jet’s lateral surface. If 
Ncr is the CR production rate in the galactic disk 
of radius Rg, the CR flux in the halo reads <I>cr — 
A’cR/(27ri?g), and the number of CRs entering the two 
jets is Nj ~ <I>cR27rA-diL?. Finally, the fraction of galac¬ 
tic CRs that can be espresso-accelerated is Nj/NcR — 
3.5%(iLj/i?g)^(Ad/2°). Extended jets with Hj ^ i?g, 
can boost a few percent of the knee nuclei by a factor 
of « 10^, producing UHECRs with energies beyond 
5Z X 10^® eV, and in particular iron nuclei with energies 
> 10^° eV. Also CR electrons, whose Galactic spectrum 
is cut off around 1 TeV, may be reprocessed via the same 
mechanism; however, it is not obvious that even TeV 
electrons have gyroradii large enough to penetrate into 
the jet and, in general, radiative losses should prevent 
them from being accelerated to very high energies. 

3.3. Spectrum and Chemical composition 

A natural prediction of the proposed mechanism is that 
the chemical composition of galactic-like CRs, which is 
increasingly heavy above 10^® eV and dominated by iron 
nuclei around 10^^eV (e.g., Kampert & Unger 2012), 
should be mapped into UHECRs. This scenario is sup¬ 
ported by recent Auger observations, which suggest a 
proton-dominated flux at 10^®eV and a heavier composi¬ 
tion at higher energies (Aab et al. 2014a). In particular, 
the composition is nitrogen-like at ^ 4x 10^® eV, whereas 
the limited statistics do not yet allow conclusive measure¬ 
ments at higher energies; the inferred trend does not ex¬ 
clude an iron contribution above 10^® eV. These results 
are not inconsistent with the Telescope Array’s report 
of a pure proton composition if the different apertures, 
event selection cuts, and Monte Carlo models are taken 
into account (Pierog 2013). With adequate statistics, the 
Telescope Array should be able to distinguish between a 
proton-only and the Auger mixed composition (see, e.g., 
Abbasi et al. 2015). 

Energy fluxes of Galactic CRs can be parameterized as 

= Ks (^^Q12 ev) (“Z«1015 ev) ’ 

where CR species (s=H, He, C/N/0, Mg/Al/Si, 
Fe) are grouped according to their (effective) charge 
Zs=l, 2, 7, 13, 26 and atomic mass As=l, 4, 
13, 27, 56; their abundances are tuned to the ones 
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E[eV] 


Figure 2. Top panel: particle fluxes above 10^® eV, for different 
species as in the legend. Dashed lines correspond to CRs acceler¬ 
ated in SNRs, and solid lines to UHECRs produced via espresso 
acceleration in AGN jets with F ^ 30. The dotted line indicates the 
extra “EeV component” (see text). Bottom panel: predicted aver¬ 
age atomic mass A as a function of energy; dashed lines correspond 
to the elements of the top panel. Colored bands represent data from 
various experiments (Kampert & Unger 2012; Gaisser et al. 2013; 
Aab et al. 2014a,b, and references therein). 


measured at 10^® eV, namely Kp ssO.lhm”^ sr“^, 
Ks/Kpi « 1,0.46,0.30,0.07,0.14, and qp Ri 2.7 and 
~ 2.6 (e.g., Horandel 2003; Caprioli et al. 2011; 
Kampert & Unger 2012). We use these simplified scal¬ 
ings, which capture the essential features of the most 
abundant species in Galactic CRs, as proxies for CR 
seeds in other galaxies as well. Interestingly, the CR 
spectrum reprocessed by AGN jets may be flatter than 
the one inside the galaxy. SNRs inject CRs in the disk 
(z = 0) with a spectrum (/>inj(E) oc E~'^, propagate in 
the halo (of thickness H Ri 3 — 5kpc) with a diffusion 
coefficient D{E) oc E^, where S r; 0.5, and escape af¬ 
ter a time Tdig ^ H^/D{E). The solution of the diffu¬ 
sion equation (e.g., Lipari 2014) returns an equilibrium 
spectrum of ()>eq oc E~'^~^ ^ E~‘^-'^ for \z\ < H, and 
(j) oc D{E)dz4'eq\z=H OC for \z\ > H. Moreover, if 
the host galaxy is very dense (such that spallation losses 
dominate over diffusive escape), one finds equilibrium 
spectra of (j)eq oc E~'^, and high-galactic-altitude spec¬ 
tra as flat as oc £’“'''+'5 (Cardillo et al. 2015). It is hence 
possible for extended jets to sweep up CR seeds with 
relatively flat spectra. Finally, energy-dependent perco¬ 
lation into the jet (and possibly seed acceleration at the 
boundary layer: see Ostrowski 1998, 2000) may lead to 
flatter injection spectra with a low-energy cut-off deter¬ 
mined by the minimum energy for which CRs can enter 
the jet. 

To calculate the spectra produced via espresso ac- 
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celeration, we take galactic-like CRs with the compo¬ 
sition in Equation 10 and with fiducial injection spec¬ 
tra (j>s{E) oc and boost them in energy by a fac¬ 

tor of Ri 10^. Note that reproducing UHECR spec¬ 
tra and composition only requires the reacceleration of 
CRs with rigidities about one decade below the knee. 
The overall UHECR normalization is chosen in order 
to reproduce the observed fluxes, and is consistent with 
the efficiency outlined above. There is a rising consen¬ 
sus (Gaisser et al. 2013; Aloisio et al. 2014; Taylor 2014) 
that an additional light component must be relevant at 
E < lO^^eV in order to explain UHECR spectra and 
composition; we consider this extra “EeV component” 
as well, modeled as in the paper by Aloisio et al. (2014). 

Figure 2 shows the obtained fluxes and chemical com¬ 
position above 10^® eV (top and bottom panels, respec¬ 
tively), compared with UHECR data. The model qual¬ 
itatively reproduces the observed spectral features and 
the characteristic light-heavy-light-heavy modulation. 
The agreement with observations is rather good, even 
if not remarkable above ^ 10^® eV. A proper description 
of source distribution, injection spectra, nuclei photo¬ 
disintegration, and proton photo-pion production would 
be needed to explain the fine structure of UHECR spec¬ 
tra and composition (Gaisser et al. 2013; Aloisio et al. 
2014; Taylor 2014), but is beyond the scope of this Let¬ 
ter. Nevertheless, we argue that inelastic losses should 
make the composition lighter at the highest energies, in 
better agreement with observations. 

4. DISCUSSION AND CONCLUSIONS 

Most of the ingredients of the presented model are 
not completely new, but we believe the proposed model 
to be quite original in its entirety. For instance, Os- 
trowski pointed out the possibility of generating UHE- 
CRs in fast blazer jets via diffusive acceleration of ther¬ 
mal particles at the jet boundary (Ostrowski 1990, 1998, 
2000). Instead, the espresso mechanism relies on ener¬ 
getic CR seeds and on a more general one-shot accelera¬ 
tion, which in addition returns the UHECR composition 
consistent with observations. Also the idea of producing 
UHECRs via CR reacceleration has already been sug¬ 
gested, either via relativistic matter plasmoids ejected 
by accreting black holes and neutron star, (“cannon¬ 
balls,” see Dar & de Riijula 2008), or via the interac¬ 
tion of AGN jets with large star-forming gaseous shells 
(e.g., Biermann et al. 2009; Gopal-Krishna et al. 2010, 
and references therein). However, such models need quite 
peculiar environmental structures whose cross-section is 
much smaller than the jet lateral surface that matters for 
espresso acceleration. 

The Lorentz factors inferred from multi-wavelength 
studies of blazars can be as high as 30-40 
(Tavecchio et al. 2010; Zhang et al. 2014). Some 
powerful sources (e.g., Mrk 421 and Mrk 501^, for 
which Lboi ^ 10"^® erg s“^ and redshift z ss 0.03) are 
found within the UHEGR “horizon,” i.e., their distance 
is less than or comparable to the mean free path for 
inelastic processes, such as photo-pion production for 
UHE protons and photo-disintegration for heavier 

^ Quite intriguingly, Mrk 421 lies within the hot-spot in the ar¬ 
rival direction of UHECRs above 57 EeV reported by the Telescope 
Array (Abbasi et al. 2014) 


nuclei (see, e.g., Dermer 2007; Aloisio et al. 2014). 
Photo-disintegration around sources can in principle 
alter the UHEGR chemical composition significantly 
(e.g. Unger et al. 2015), but it is not expected to be 
important for the espresso mechanism, if most of the 
reacceleration occurs well beyond the AGN broad line 
region (see, e.g., Dermer 2007). 

The low level of anisotropy measured in the UHEGR 
arrival directions (Abreu et al. 2010; Abbasi et al. 2014) 
suggests that the horizon encompasses several sources, 
not necessarily associated with copious emission in other 
bands, such as 7 —rays. The census of potentially- 
relevant active galaxies may not be limited to the few 
powerful known blazars inside the horizon; since some 
espresso-accelerated UHECRs may escape the jet side¬ 
ways, also radio galaxies —whose 7 -ray emission is not 
beamed in our direction— may partially contribute to 
the measured UHECR flux. 

In general, different types of AGNs (not just powerful 
BL Lac objects and flat-spectrum radio quasars) may be 
responsible for the global shape of the galactic to extra- 
galactic transition. The EeV component, for instance, 
may be produced via espresso acceleration in galaxies 
hosting jets with moderate luminosities and T ;< 5 (see 
Equation 1), which are numerous inside the horizon. 
Such AGNs could be responsible for the production of the 
multi-PeV neutrinos detected by IceCube (Aartsen et al. 
2013). A more detailed account for AGN zoology (dis¬ 
tribution in luminosity, Lorentz factors, and redshift) is 
needed to further probe this scenario. 

Finally, we check whether espresso acceleration may 
work also in the ultra-relativistic flows of other candidate 
sources of UHECRs. Pulsar winds have very large P > 
10®, and live in environments rich in CR seeds accelerated 
at the surrounding SNR shocks. However, for the Crab 
parameters (P ss 10®, D « O.lpc, and B' ^ Bpwn/P, 
where HpwN ~ 100/iG is the nebular field), CRs cannot 
achieve isotropization within the wind (see Equation 7); 
such a requirement is unlikely to be fulfilled also in newly 
born millisecond pulsars. For GRBs, instead, the limited 
extent of their jets makes it impossible to reaccelerate a 
sizable fraction of galactic CR seeds. 

In summary, we presented a general mechanism for the 
acceleration of UHECRs in energetic AGN jets. If the 
relativistic flow is sufficiently extended (Equations 2 and 
7), seed CRs with gyroradii large enough to cross the jet 
boundary can be boosted by a factor of ~ P^ in energy 
and be accelerated up to 10^° eV. Under reasonable 
assumptions, such a one-shot espresso acceleration can 
account for both the spectrum and chemical composi¬ 
tion of UHECRs (Figure 2). In forthcoming works we 
will probe this scenario by: testing CR injection and en¬ 
ergy gain in realistic jet structures; accounting for AGN 
phenomenology; and including UHEGR propagation and 
losses. 
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